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Abstract 



Altruistic behaviour is disadvantageous for the individual while is advanta- 
geous for its group. If the target of the selection is the individual, one would 
expect the selection process to lead to populations formed by wholly homo- 
geneous groups, made up of either altruistic or egoistic individuals, where the 
winning choice depends on the balance beetwen group advantage and indi- 
vidual disadvantage. We show in a simple model that populations formed by 
inhomogeneous groups can be stabilized in some circumstances. We argue 
that this condition is realized when there is a relative advantage conferred by 
the presence of a few altruists to all the members of the group. 
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I. INTRODUCTION 

Natural selection is an egoistic process: individual organisms compete for representation 
of their genes in the next generation, and only those genomes better able to survive and 
reproduce are likely to be maintained in the evolutionary process. In this framework it 
is difficult to explain how "altruistic" genes (that is genes which determine a behaviour 
disadvantageous to the carrier, but beneficial for other individuals) can be selected. 

Nevertheless, we can find in nature some examples of altruistic behaviour: social or- 
ganization in insect societies H0, parental care |]9|JTT[1, warning calls in birds fllQfl , and so 



on. So far, there is no evidence of a genetical basis for altruistic behaviour. However, it is 
an interesting problem in itself to understand how interactions between individuals may be 
reflected on the action of individual natural selection. The presence of altruists increases 
the chances of the group which it belongs to. The general problem is not only to understand 
what is the level of selection action, but also how the different levels may interact under the 
selective pressure. 

Classically, two main mechanisms have been proposed to explain the origin and evolution 
of altruism: kin selection [|], ||[7]], [PkPH, and group selection fl^fl , [[|. In both mechanisms 



it is necessary to assume that individual fitness depends not only on individual genotype, 
but also on the genetical composition of the population. In kin selection the individual 
probability of survival is proportional to the relatedness that the individual shares with 
the altruist. The selection mechanism thus is limited only to related individuals, and the 
selection unit is the individual. In group selection, the interaction takes place between 
individuals belonging to the same group, and groups are defined by external constraints 
(for example spatial boundaries), disregarding any relatedness between individuals. In most 
cases, one ends up with a situation in which either egoistic or altruistic behavior is selected. 
Nevertheless, it is clear from observation that most often altruistic behavior is exhibited 
only by some individuals in the group. 

A recent report concerns territorial defense by prides of lions. Female lions living in 
Serengeti National Park and Ngorongoro Crater (Tanzania) live in groups that hunt togheter 
and defend a common territory against other groups. Heinsohn and Packer report in || 
that some lionesses act more aggressively and swiftly in defense, while others tend to lag 
behind. No consistent relations were observed between physical size and defense-readiness. 
Furthermore, when engaging in aid of a pride-mate mounting a defense, no consistent relation 
was found between the readiness of the helper, and the kin relation between the helper and 
the first defender. Territorial defense in groups of female lions therefore seems to be a case 
of altruistic behaviour where kin relation does not play a role. 

We show in the following, in an exactly solvable model, that altruistic behavior can be 
selected both disregarding relatedness between individuals and maintaining the individual as 
the selection unit: morevoer, that, if the advantage of the presence of altruistically behaving 
individuals is already felt when their number is small, the stable situation is one in which 
inhomogeneous groups, formed of both egoistic and altruistic individuals, coexist with both 
fully egoistic and fully altruistic individuals. 
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TABLE I. Fitness table. Rows: individual genotypes; columns: genetical group composition. 



We consider a population made up of M individuals, divided in N groups (denies) of 
L = M/N individuals on average. We assume non-overlapping generations, heredity acting 
according to the usual Mendelian mechanism, and the presence of a behavioral locus with 
two alleles: A (recessive) and E (dominant). An individual of genotype AA is altruist, while 
AE and EE are not. However, as a first approximation, we disregard the individual genetical 
structure, considering only the two classes of individuals: altruist (A) and egoist (E). Demes 
have a variable number j of altruistic individuals A. If, for a given deme, j is greater than 
a threshold j* the deme is called altruistic, otherwise egoistic. We shall call any deme with 
exactly j altruistic individuals a j-deme. 

To each individual a in a generation we associate its fitness value f(a), which is propor- 
tional to the survival probability of its genome in the population over one step. The fitness 
f(a) of the individual a depends on two factors: 

1. on its genotype: if the individual is altruist (genotype A), its fitness is reduced by a 
factor (1 — r) with respect to the other members of its deme; 

2. on its deme composition: if the individual (whatever its genotype) belongs to an 
altruistic deme its fitness is enhanced by a factor 1/(1 — c), with respect to individuals 
with the same genotype belonging to an egoistic one. 

The parameter r is called intrademic selection rate and c is called interdemic selection rate. 
This form of the fitness is summarized in Table 1. 

The population evolves according to a discrete process, involving the following steps: 

1. Selection and reproduction: Taking into account the internal genetical structure, one 
can consider both asexual and sexual reproduction mechanisms. In the following, we 
will consider only the asexual case, but the reproduction mechanism does not affect 
the essential results. The average number of the offspring of an individual a is given 



where f(a) is the fitness of the individual a, as defined in Table [IJ and the sum runs 
over the whole population. At the end of this step, the offspring of the individuals of 
one deme in one generation form one deme. 

2. Deme splitting: The mean size of the demes is equal to L = M/N . However, after 
the reproduction step, some demes can disappear and the mean size of the surviving 



II. THE MODEL 
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demes would correspondingly increase. In particular, some denies may have become 
much greater than L. We assume that too large demes split in such a way as to keep 
the mean size L constant. This rule is motivated by the analogy with social animals, 
which do not typically live in too large groups. It is also a necessary ingredient for 
purely modelling reasons: if there is no maximal size of demes, sooner or later one of 
the demes would take over the whole population, and interdemic selection would cease 
to operate. 

Notice that an A individual has always less chances than an E individual belonging to 
the same deme, but all the individuals in an altruistic deme may have a greater survival 
probability than individuals belonging to egoistic demes, independently of their genotypes. 



III. THE MASTER EQUATION 

We assume that the total number of individuals M tends to infinity along with deme 
number N, keeping the typical deme size L finite. We denote by nj(t) the number of j-demes 
at generation t, and by Pj(t) = rij(t)/N the corresponding fraction. The fractions Pj{t) are 
normalized: 

Eft-W =1, v*- (2) 

j 

Equation ([!]) yields the average number of offspring of a single individual. Assuming that 
reproduction events are independent from one another, the number of offspring of the single 
individual is Poisson distributed with average equal to the expected value v{pt) (eq. (0)). As 
a consequence, the number of offspring of a whole deme S is also Poisson distributed, with 
average Q(S) given by the sum of v(a>) over its members: 



Q05) = 5>(a). (3) 

It follows from Eq. ([J) that the size of the offspring of the j'-deme 5 depends only on the 
number j of its altruists and is given, on average, by Q(S) = VjL, where 

V ^ j 4 

In this equation we have introduced the function 

_ 1 ~ rj/L 

3 ~ l-c9(j-j*) ' {b) 

where 6{x) is the unit step function. We stress that this is the average new size of the deme, 
disregarding its new genetic composition. Indeed, the genetic composition of a new deme 
which originates from a j-deme is also random. We recall that the number of altruists in 
the old deme is j while the corresponding number of egoists is L — j. We also recall that 
the relative probability of survival in the same deme is 1 — r. The probability that a given 
offspring of a j-deme is altruist is therefore 
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Pi = -J— ■ ( 6 ) 

Since reproduction events are independent, the number of altruists is binomially distributed 
with success probability p = Pj. This means that the average fraction of altruists is equal 
to pj, while its variance is approximately given by 

T P = Pi(l-p i )^ E . (7) 

This expression holds in the approximation in which the average size of the offspring deme 
is substituted for its actual size: due to the reproduction, demes grow up and shrink, so 
that their actual size is not constant. 

We now consider the effects of deme splitting. It produces new demes with smaller 
size but with the same composition in average. Assuming independence, we have again a 
binomial redistribution of individuals in the splitted demes. The average size is brought 
back to L and the average fraction of altruists is still given by pj, while its variance is given 
by 



a 



split - '■ 



J 



P J (1-P J ) I . (8) 



If the deme size L is large enough, we can approximate the distribution of altruist 
numbers in the new demes by a gaussian. Define 



where 



<7 i = O f Ut + 0? P . (10) 



The transition probability matrix (i.e., the probability of having i altruist in a deme at time 
t + 1 from a j-deme at time t, after the reproduction and splitting steps) is given by 

A ii = ^- i (j^0,L). (11) 

Since, in the absence of mutations, all demes originating from a homogeneous deme (be it 
altruistic or egoistic) are also homogeneous and of the same type, we have 

(12) 
(13) 
(14) 

It should be remarked that only large demes split and therefore <jj = erj ep for small demes. 
Nevertheless the error we introduce extending (|T0D to all demes is negligible, because for 
small demes Vj < 1, so that cr| pllt < o\j cp . 



A LL 


= A 00 


= 1; 


A iL 


= 0, 




A i0 


= 0, 


i ^ 



5 



The number of denies changes because of splitting: however, since the number of demes, 
N, goes to infinity, the fraction of demes originating from j-demes is equal to its average 
because of the law of large numbers. This number is given by 

v '* w = f£&' (15) 

We are thus led to the equation 

ft^+^^EM-ftW » ( 16 ) 

which closely resembles Eigen's quasispecies equation 0. The normalization factor Z t is 
given by 

Z(t) = I2T jPj (t) (17) 

j 

The main difference between eq. ([16]) and the quasispecies equation lies in the fact that the 
matrix of transition probabilities (Aij) also depends on the vector (pj(t)) via the variances 

to)- 



IV. RESULTS 

In the long run, the population will reach a steady distribution, which may be of one 
among three types. The first two types correspond to a population composed only by 
altruistic or egoistic individuals. It is trivial to remark that the demes for these two phases 
have an homogeneous composition. The third type correspond to a situation where, together 
with homogeneous demes of both kinds, one has inhomogeneous demes, formed by both 
egoistic and altruistic individuals. 

Let us first imagine that the population is initially composed only by homogeneous demes 
of either kind: completely altruistic and completely egoistic. In this case no mixed demes 
can be generated, since homogeneous demes only reproduce themselves. The advantage for 
the completely altruistic deme is 1/(1 — c) while its disadvantage is 1 — r. It is easy to 
understand that, in this case, the population will end up being made of only completely 
altruistic demes if c > r and only egoistic ones if c < r. In other words: 

^' )= {J; foIi>o| forc<r; (i8a) 

Pe q (j) = { tVj^L; fOTC>r - (18b) 

If we mark with dots the region in which only altruistic individuals are present and with short 
vertical lines the region where only egoistic ones are present in the (r, c) plane, we obtain 
the "phase diagram" of Fig. 1. We call "phases" the different regimes of the equilibrium 
population: egoistic, altruistic, or mixed. 
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FIG. 1. Phase diagram in the (r, c) space for L (deme size) equal to 20 and j* (threshold) 
equal to 4. One observes the stability of a wholly altruistic population (dots) for c > r and a fully 
egoistic one (vertical lines) for c < r. 

Suppose now that also inhomogeneous denies are initially present. In this case that 
denies with j > j* share with completely altruistic ones the "altruistic advantage" and, 
since they have L — j egoistic individuals, they partially share the "egoistic advantage" with 
the completely egoistic demes. From this point of view, the fittest demes would be those 
with j = j*. Nevertheless, reproduction and splitting changes their composition, so that 
they most of the times do not reproduce. Briefly, mixed demes can generate homogeneous 
demes but not viceversa. Due to these competing effects, it not a priori clear if the dynamics 
may lead to mixed populations. If it is so, we expect that this will happen provided that j* 
is so small that the mixed demes can take advantage of both altruistic and egoistic behavior 
in the best way. This is shown by the equilibrium solution of the master equation. Compare, 
in fact, Fig. 1 with Fig. 2, in which L = 20 again, but j* = 1 instead of 4. One sees the 
co-existence, at equilibrium, of both altruistic and egoistic individuals. 

Figs. 3 and 4 correspond to intermediate situations with j* = 2 and j* = 3. It is clear that 
the region where mixed populations are stable becomes progressively smaller from j* = 1 to 
j* = 3 and disappears for j* = 4. 

Let us stress that the equilibrium distribution p eq (j) is non-trivial when the equilibrium 
population is inhomogeneous (see, for example, Fig. [|). In fact, together with mixed demes, 
also homogeneous demes are present since they arise as offsprings of mixed demes at any 
generation. The edges of the stability region of the mixed phase depend both on L and on 
j*. Figures 1-4 correspond to L — 20, but for larger L the qualitative behaviour it is the 
same: if we increase j* the domain of mixed phase becomes smaller and, above a critical 
value (jc r ), it desappears. The value of j* r depends on L. Plotting j* r vs. L, we obtain the 
graph represented in Fig. [|, from which one clearly sees that 

hm j c ;/L =0, (19) 

L— >oo 

while j* r /L = const., for L < 60. This behaviour suggests that for large demes the mixed 
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FIG. 2. Phase diagram in the (r, c) space for L = 20 and j* = 1. One observes the appearance 
of a region of stability of an inhomogeneous (mixed altruistic-egoistic) population (diamonds). 
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FIG. 3. Phase diagram in the (r, c) space for L = 20 and j* = 2. The region of stability of the 
inhomogeneous population (diamonds) has shrinked. 
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FIG. 4. Distribution of the number j of individuals in each deme for L = 20, j* = 2, r = 0.4, 
c = 0.0, 0.05, 0.1, . . . , 0.95. One sees that most demes are inhomogeneous (j ^ 0, L) and none is 
formed only by altruistic individuals. 

phase can exist only if the altruistic advantage is granted by a small number of altruistic 
individuals. 

V. DISCUSSION 

The most important feature of our model is that selection acts at individual level, al- 
though the individual survival probability takes into account also of the "interaction" among 
individuals. This very simple interaction is modeled by the dependence of the fitness func- 
tion (see Tab. 1) on the genetic composition of group: the individual survival probability 
is sensitive to the presence of other altruists in the group. So, without invoking either a 
kin selection mechanism or a group selection one for the maintenance of altruists in the 
population, we have seen that, when the group advantage balances the individual (altruis- 
tic) disadvantage, mixed groups may co-exist at equilibrium. This kind of equilibrium also 
depends on the percentage of altruists in a deme (for small enough demes, i.e., L < 60); for 
larger demes a finite number of altruists is sufficient to produce inhomogeneous demes at 
equilibrium. 

Let us remark that the stabilization of inhomogeneous populations is reached in the 
absence of the (rather improbable) mutations of the "altruistic" into the "egoistic" gene 
(and viceversa): the key role is played by the fluctuations in the composition of daughter 
demes from inhomogeneous ones. These fluctuations are even more important when, as in 
Ref. [IIJ, the size of the daughter demes is allowed to vary. 

Summarizing, we have shown in a numerically solvable model that if the advantage of 
the altruistic behavior of a few individuals falls (above a given threshold) on the whole of a 
sufficiently large deme, a stable situation arises in which inhomogeneous demes (consisting of 
both altruistic and egoistic individuals) coexist with homeogeneous ones. For large average 
deme size, the average number of altruistic individuals approaches the threshold. 
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FIG. 5. The critical value j* r of j* , for which the population becomes homogeneous, as a 
function of deme size L. For L < 60, one has j* r oc L, while j* r attains a constant value for L > 60. 
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